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We study electron transport in quasi-one-dimensional metallic wires. Our aim is to compare an impurity- 
free wire with rough edges with a smooth wire with impurity disorder. We calculate the electron transmission 
through the wires by the scattering-matrix method, and we find the Landauer conductance for a large ensemble 
of disordered wires. We first study the impurity-free wire whose edges have roughness with a correlation length 
comparable with the Fermi wave length. Simulating wires with the number of the conducting channels (N c ) 
as large as 34 - 347, we observe the roughness-mediated effects which are not observable for small N c (~ 3 
- 9) used in previous works. First, we observe the crossover from the quasi-ballistic transport to the diffusive 
one, where the ratio of the quasi-ballistic resistivity to the diffusive resistivity is ~ N c independently on the 
parameters of roughness. Second, we find that transport in the diffusive regime is carried by a small effective 
number of open channels, equal to ~ 6. This number is universal - independent on N c and on the parameters 
of roughness. Third, we see that the inverse mean conductance rises linearly with the wire length (a sign of 
the diffusive regime) up to the length twice larger than the electron localization length. We develop a theory 
based on the weak-scattering limit and semiclassical Boltzmann equation, and we explain the first and second 
observations analytically. For impurity disorder we find a standard diffusive behavior. Finally, we derive from 
the Boltzmann equation the semiclassical electron mean-free path and we compare it with the quantum mean- 
free path obtained from the Landauer conductance. They coincide for the impurity disorder, however, for the 
edge roughness they strongly differ, i.e., the diffusive transport in the wire with rough edges is not semiclassical. 
It becomes semiclassical only for roughness with large correlation length. The conductance then behaves like 
the conductance of the wire with impurities, also showing the conductance fluctuations of the same size. 

PACS numbers: 73.23.-b, 73.20.Fz 



I. INTRODUCTION 



A wire made of the normal metal is called mesoscopic if the 
wire length (L) is smaller than the electron coherence length 
It is called quasi-one-dimensional (QID), if L is much 
larger than the width (W) and thickness (H) of the wire [3]. 
Fabrication of the QID wires from such metals like Au, Ag, 
Cu, etc., usually involves techniques like the electron beam 
lithography, lift-off, and metal evaporation. These techniques 
always provide wires with disorder due to the grain bound- 
aries, impurity atoms and rough wire edges J3l. Disorder scat- 
ters the conduction electrons and limits the electron mean free 
path (I) in the wires to ~ 10 — lOOnm Of fundamental in- 
terest are the wires with W and H as small as ~ 10 — lOOnm. 

In this work the electron transport in metallic QID wires 
is studied theoretically. We compare an impurity-free wire 
with rough edges with a smooth wire with impurity disorder 
(a wire with grain boundaries will be studied elsewhere). We 
study the QID wires made of a two-dimensional (2D) conduc- 
tor (H — > 0) of width W and length L > W. Our results are 
representative for wires made of a normal metal as well as of 
a 2D electron gas at a semiconductor heterointerface. 

First we review the basic properties of the QID wires. Con- 
sider the electron gas confined in the two-dimensional (2D) 
conductor depicted in the figure [T] At zero temperature, the 
wave function (p(x, y) of the electron at the Fermi level {Ep) 
is described by the Schrodinger equation 




FIG. 1: Wire made of the 2D conductor of width W and length L. 
The figure on the left depicts the wire with impurities positioned at 
random with random signs of the impurity potentials. The figure on 
the right depicts the impurity-free wire with rough edges, where d(x) 
and h(x) are the y-coordinates of the edges at y = and y = W, 
respectively, randomly fluctuating with x. The roughness amplitude 
is A and the step Aa; also means the correlation length (see the text). 



with Hamiltonian 

where m is the electron effective mass, Ui(x, y) is the poten- 
tial due to the impurities, and V(x,y) is the confining poten- 
tial due to the edges. Following the figure Q] the confining 
potential in a wire with smooth edges can be written as 



V(y) 



0, < y < W 
oo, elsewhere 



while in a wire with rough edges it has to be modified as 



Hip{x,y) = E F tp(x,y) 



(1) 



V(x,y) 



0, d(x) < y < h(x) 
oo, elsewhere 



(3) 



(4) 
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where d(x) and h(x) are the y-coordinates of the edges. The 
potential of the impurities, Uj, is usually assumed to be a 
white-noise potential [3]. The simplest specific choice is 



(5) 



where one sums over the random impurity positions [xi,yi] 
with a random sign of the impurity strength 7 (see Fig. [JJ). 
Similar models of disorder like in the figure [JJ are commonly 
used in the quantum transport simulations Bq- tioll 

The disordered QID wire is connected to two ballistic semi- 
infinite contacts of constant width W, as shown in the figure 
12 In the contacts the electrons obey the Schrodinger equation 



"2m \dx 2 + dy 2 * (y > 



(p{x,y) = E F ip(x,y), 



(6) 

where V(y) is the confining potential given by equation ([3j. 
Solving equation (O one finds the independent solutions 



<pj;(x,y)=e ±ik »*xM n=l,2, 
with the wave vectors k n given by equation 

H 2 kl h 2 n 2 



. . 00, 



Ep 



2m ' 2mW 2 ' 

where e n is the energy of motion in the j/-direction and 



Xn(y) 



0. 



y), < y < W 
elsewhere 



(7) 



(8) 



(9) 



is the wave function in the direction y. The vectors k n in (JTJ 
are assumed to be positive, i.e., the waves e tknX and e~ lk " x 
describe the free motion in the positive and negative direc- 
tion of the x-axis, respectively. The energy e„ + h 2 k 2 /2m is 
called the n-th energy channel. The channels with e n < Ep 
are conducting due to the real values of k n while the chan- 
nels with e n > Ep are evanescent due to the imaginary k n . 
The conducting state <p^(x, y) = e lknX Xn(y) in the contact 1 
impinges the disordered wire from the left. It is partly trans- 
mitted through disorder and enters the contact 2 in the form 



<P%(x,y) = tmn e * kmX x m (y), x > L, 



(10) 



where t mn (k n ) is the probability amplitude of transmission 
from n to m. At zero temperature, the conductance of the 
disordered wire is given by the Landauer formula 111 ] 



G 



2e 2 ^ 

n=l 



__ \ " \ " U |2 ^ 



n—l m—1 



kn 



(ID 



where we sum over all (N c ) conducting channels. We note 
that T n is the transmission probability through disorder for 
the electron impinging disorder in the n-th conducting chan- 
nel. The amplitudes t rnn have to be calculated for specific dis- 
order by solving the equation ((TJ with the asymptotic condi- 
tion ([Tol l. In the ensemble of macroscopically identical wires 



contact 1 



contact 2 





QID wire 


IF 
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FIG. 2: The QID wire placed between two contacts. The bold 
arrows denote the wave amplitudes A + , B coming in the wire 
and the amplitudes A ~ , B + coming out the wire. 



disorder fluctuates from wire to wire and so does the conduc- 
tance. Hence it is meaningful to evaluate (fTTb for the ensem- 
ble of wires and to study the ensemble-averaged conductance 
(G), variance (G 2 ) — (G) 2 , resistance (1/G), etc. We now 
discuss a few important results of such studies. For simplicity 
we use the variables g = Gj (2e 2 /h) and p = 1/g. 

First we discuss the smooth QID wires with disorder due 
to the white-noise potential Ui. For L — the formula (fTTb 
gives the ballistic conductance g = N c . As L increases, the 
formula ( fTTT i first shows the classical transmission law [2] 



(.9) = ^Vc 



-I 



< L<£, 



(12) 



where N c ^S> 1 and £ ~ NJ is the QID localization length. If 
I <C L <C £, the wire is in the diffusive regime. For / <C L and 
N c ~ kpWI-K we obtain from ( TT2l the standard expression 

id) =<7difW/L, a dif = nn e l/k F , I < L < £, (13) 

where o^ij is the diffusive conductivity, kp is the 2D Fermi 
wave vector, and n e = k 2 F /2ir is the 2D electron density. 
However, the mesoscopic diffusive conductance is also af- 
fected by weak localization. Hence, one in fact obtains from 
(fTTl i a slightly modified version of (fT3l , namely [J, 1 1 2TI 



(9) =a dif W/L- 1/3, ?«L<C£, 



(14) 



where the term 1 /3 is the weak localization correction typical 
of the QID wire. The mean free path I in the above formulae 
coincides with the mean free path derived from the semiclas- 
sical Boltzmann transport equation, i.e., the quantum conduc- 
tance (fTTT i captures the Boltzmann transport limit exactly. The 
mean of the two-terminal resistance, (p) = (1/g), shows in 
absence of weak localization the diffusive behavior yO 



(p) = p c + p d ifL/W, 0<L«£, 



(15) 



where p c — 1/N C is the contact resistance and pdif = 1/vdif 
is the diffusive resistivity. The diffusive resistance dTsb and 
diffusive conductance ( fT3b thus coexist in a standard way: 
(p) ~ I /(g) for Pdifj^r Pc- If we include the weak lo- 
calization by means of ( TT4"1 >. then 



(P) - l/(.9} - Pdlfyy + 



1 2 I?_ 



/<L<£. (16) 
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The conductance fluctuates in the diffusive regime as fl 1 3L 11411 

Vvar( 5 ) = V(.9 2 ) - <5> 2 = VVW ~ 0.365, k< L « ^, 

(17) 

where the numerical factor 0.365 is typical of the Q1D wire. 
Finally, as L exceeds £, the mesoscopic Q1D wire enters the 
regime of strong localization, where (g) decreases with L ex- 
ponentially while (p) shows exponential increase IU5U16I1 . 

The formulae ( flZb - ST% hold for the wires with impurity 
disorder. Do they hold also for the wires with rough edges? 
In our present work we address this question from the first 
principles: we calculate the amplitudes t mn by the scattering- 
matrix method 0, for a large ensemble of macroscop- 
ically identical disordered wires, we evaluate the Landauer 
conductance (fTTt . and we perform ensemble averaging. 

In fact, a few serious differences between the wires with 
rough edges and wires with impurity disorder were identified 
prior to our work. In the wire with impurity disorder the chan- 
nels are equivalent in the sense that (Ti) = (T2) • • • = (Tjv c ) 
1 18, 19ll . For instance, in the diffusive regime (T n ) = \l/L 
for all channels |2|]. In the impurity-free wires with rough 
edges (T n ) decays fast with rasing n, because the scattering by 
the rough edges is weakest in the channel n = 1 and strongest 
in the channel n = N c Il20l - l24ll . This is easy to understand 
classically: in the channel n = 1 the electron avoids the edges 
by moving in parallel with them, while in the channel n = N c 
the motion is almost perpendicular to the edges, resulting in 
frequent collisions with them. As a result, (T n ) shows coex- 
istence of the quasi-ballistic, diffusive, and strongly-localized 
channels II25I1 . Due to this coexistence, the work 12011 reported 
absence of the dependence (g) oc 1/L, suggesting that the 
wire with rough edges does not exhibit the diffusive conduc- 
tance ( fT3l . However, according to Il2ll 12611 . the wire with 
rough edges seems to exhibit the diffusive resistance ( fT~5b . In 
our present work these findings are examined again, but for 
significantly larger N c as in previous works. 

First we study the impurity-free wire whose edges have a 
roughness correlation length comparable with the Fermi wave 
length. For L — > we observe the quasi-ballistic dependence 
l/{g) = (p) = l/N c + p q bL/W, where p q b is the quasi- 
ballistic resistivity. As L increases, we observe crossover to 
the diffusive dependence l/(g) ~ (p) = l/7V c e// +p dlf L/W, 
where pdi / <C p q b and 1 /N*** is the effective contact resis- 
tance due to the N*** open channels. We find the universal 
results pqb/ Pd t f ^ 0.6iV c and N^ ff ~ 6 for N c > 1. As L 
exceeds the localization length £, the resistance shows onset 
of localization while the conductance shows the diffusive de- 
pendence 1/(5) ~ l/iV c e// + pdif L/W up to L ~ 2£ and the 
localization for L > 2£ only. Finally, we find 



v/var(.g) = V V> - (g) 2 ^ 0.3, I <C L < 2£. (18) 

The fluctuations ( fT8l differ from (TTTb and were already re- 
ported in the past Q2L 1^51 127L [2^] . For the smooth wires with 
impurities our calculations confirm the formulae (fT2l - dTTb . 

Moreover, we derive the wire conductivity from the semi- 
classical Boltzmann equation 029143111 . and we compare the 



semiclassical mean-free path with the mean-free path ob- 
tained from the quantum resistivity pdif- For the impurity 
disorder we find that the semiclassical and quantum mean- 
free paths coincide, which is a standard result. However, for 
the edge roughness the semiclassical mean-free path strongly 
differs from the quantum one, showing that the diffusive trans- 
port in the wire with rough edges is not semiclassical. We 
show that it is semiclassical only if the roughness-correlation 
length is much larger than the Fermi wave length. For such 
edge roughness the conductance behaves like the conductance 
of the wire with impurities (formulae [T2l - \T%, also showing 
the fluctuations (ITTb . 

The next section describes the scattering-matrix calcula- 
tion of the amplitudes t mn for the impurity disorder and edge 
roughness. In section III, the impurity disorder and edge 
roughness are treated by means of the Boltzmann equation 
and the semiclassical Q1D conductivity expressions are de- 
rived. In section IV we show our numerical results. More- 
over, the crossover from I /(g) = (p) = 1/N C + p q bL/W to 
I /(g) ~ (p) = 1/A r c e// + pdif L/W in the wire with rough 
edges is derived by means of a microscopic analytical theory. 
The theory neglects localization but nevertheless captures the 

main features of our numerical results. In particular, we ob- 
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tain analytically the universal results p q b/ Pdif = fj^c and 
Ngff ~ 2.5. A summary is given in section V. 



II. THE SCATTERING-MATRIX APPROACH 

Consider the Q1D wire with contacts 1 and 2, shown in the 
figure |2 The wave function ip(x,y) in the contacts can be 
expanded in the basis of the eigenstates ©. We introduce no- 
tations A^(x) = a^e ±tk " x and -B r f (x) = b^e ±lk " x , where 
and are the amplitudes of the waves moving in the pos- 
itive and negative directions of the x axis, respectively. At the 
boundary x = 

N 

<p(0,y) = ]T [A+(0) +A-(0)] X n(y), (19) 

n=l 

while at the boundary x = L 

N 

<p(L,y) = V [B+(L)+B-(L)] X n(y), (20) 



where N is the considered number of channels (ideally N — 
00). We define the vectors A ± (0) and B ± (L) with compo- 
nents Ar_, jy(O) and -B„ =1 n(L), respectively, and we 
simplify the notations A ± (0) andB ± (i) as A" 11 andB"* 1 . The 
amplitudes A 1 * 1 andB^ 1 are related through the matrix equation 



(21) 



where S is the scattering matrix. Its dimensions are 2A^ x 2A^ 
and its elements t, r, t', and r' are the matrices with dimen- 
sions N x N, Physically, t and t' are the transmission am- 
plitudes of the waves A + and B , respectively, while r and 





~ r t'~ 




, S = 


' r f ' 


(SO- 


t r' 


(»-) 


t r' 
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r' are the corresponding reflection amplitudes. The matrix el- 
ements of the transmission matrix t are just the transmission 
amplitudes t mn which determine the conductance dTTb . 

Consider two wires 1 and 2, described by the scattering ma- 
trices S\ and S 2 . The matrices are defined as 



Si = 



Let 



n A 



S 2 = 



r 2 t> 2 
h r' 



ri2 t\ 



12 



(22) 



(23) 



is the scattering matrix of the wire obtained by connecting 
the wires 1 and 2 in series. The matrix S12 is related to the 
matrices Si and S 2 through the matrix equations [2] 



tia =t 2 [I-r' 1 r 2 ]-H 1 , 
ri2 =n+ t[r 2 [I - r' l r 2 ]~ 1 tii 
t' 12 =t' 1 [I + r 2 [I-r[r 2 ]- l r[]t' 2 , 
r' 12 =r' 2 + t 2 [I-r' 1 r 2 ]- l r' l t l 2 ,, 



(24) 



' 2" 



where / is the unit matrix. The equations d24b are usually 
written in the symbolic form 



Si 2 = Si®S 2 . 



(25) 



A. Scattering matrix of smooth wire with impurity disorder 

Consider the wire with impurity potential ©. Between any 
two neighboring impurities there is a region with zero impu- 
rity potential, say the region Xi-i < x < Xi, where the elec- 
tron moves along the x axis like a free particle. The wire with 
n impurities contains n + 1 regions with free electron motion, 
separated by n point-like regions where the scattering takes 
place. As illustrated in figure [3] the scattering matrix S of 
such wire can be obtained by applying the combination law 



S = pi (g> Si ®p 2 ® S 2 



'Pn+l, 



(26) 



where pi is the scattering matrix of free motion in the region 
x^i < x < Xi and Si is the scattering matrix of the i-th 
impurity. The symbols eg) mean that the composition law d2~5T l 
is applied in ( f26b step by step: one first combines the matrices 
pi and si, the resulting matrix is combined with p 2 , etc. 
The scattering matrix pi can be expressed as 



Pi 



$ 
$ 



(27) 



where is the N x N matrix with zero matrix elements and 
$ is the N x N matrix with matrix elements 

^mn — ^ &mn: Ci — 1; (28) 

Finally, for a <5-function-like impurity the scattering matrix 

r t' ' 

t r> ^ 



1 


c c 

m 1 ^ a 2 


2 




+4 w 




A + 


\ -4 w 


B + 






A 




B 



y\ , p 1 ® s ® p 2 9 s 2 

X 



FIG. 3: Wire with the randomly positioned point-like impurities. The 
n impurities described by the scattering matrices s» divide the wire 
into the n + 1 free regions described by the matrices pi. Also shown 
are the wave amplitudes A ± andB*. 



is composed of the matrices 10,0] 



t = t' = [K + iT^K, 
r = r' = -[K + iT]- 1 iT, 



(30) 
(31) 



where K and Y are the N x N matrices with matrix elements 

K mn = k n S mn , T mn — "^2"X m (j/i)Xn(j/i)- (32) 

Concerning the value of N, we use N > N c chosen in such 
way that in the diffusive regime our simulation reproduces 
the Boltzmann-equation results. 



B. Scattering matrix of the impurity-free wire with rough edges 

The electrons in the impurity-free wire with rough edges are 
described by the Schrodinger equation (Q]) with Hamiltonian 
without the impurity potential, but with the confining poten- 
tial V(x, y) [equation ©] including the edge roughness. We 
specify the edge roughness as follows. We define Xj = jAx, 
where j — 0,1,2,... and Ax is a constant step. For x be- 
tween xj and Xj+i, the smoothly varying functions V(x, y), 
h(x), and d(x) in equation are replaced by constant values 
Vj(y) = V(xj,y), hj = h(xj), and dj = d(xj), respectively. 
We obtain the equation 



Vj(y) 



o, 

00, 



dj < y < hj 
elsewhere. 



(33) 



We assume that dj and hj vary with varying j at random in 
the intervals (—A, A) and (W — A,W + A), respectively. 
This is depicted in the figureQ] where h(x) and d(x) fluctuate 
with varying x by changing abruptly after each step Ax. 
The wire width fluctuates with varying x as well. However, 



for x between Xj and Xj+i we have the constant width 



Wi 



(34) 
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p®s®p®S®... ®o 



n+] 



FIG. 4: The impurity-free wire with n edge steps described by the 
scattering matrices Sj and with n + 1 free regions of constant width 
Wj, described by the scattering matrices Pj. 
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FIG. 5: The wire with a single edge step at x = 0. The symbols in 
the figure are discussed in the text. 



Consequently, the electron wave function for x between Xj 
and Xj + i can be expressed in the form 



JV, 



(x,y) = J2[a+e ik > 



a„ e 



-ik-L x 



xUv), (35) 



where Nj is the considered number of channels (Nj = oo in 
the ideal case), the wave vectors k 3 n are given by equations 



2m 



xUv) 



(y - dj) 



o. 



h 2 ir 2 2 



dj < y < hj 
elsewhere 



(36) 



(37) 



are the wave functions for the y-direction, and the index j 
means that the above equations hold for x between Xj and 
Xj+i. The equations d36t and d37l i are just the equations dHJ 
and (0, respectively, modified for the wire with rough edges. 
In practice we choose Nj by means of the relation lfl7ll 



N j /N = Wj/W, 



(38) 



i.e., the ratio of the channel numbers in the j-th region and 
contact regions is the same as the ratio of their widths. We 
choose N > N c and we check, that the calculated conduc- 
tance does not depend on the choice of N. For N > N c the 
relation d38l > ensures N^ > Nj!, where Nj! is the number of 
the conducting channels in the j-th wire region. This makes 
the calculation reliable also when N[ fluctuates with varying 
j, which happens for A larger than the Fermi wave length. 

As shown in figure @1 the scattering matrix S of the wire 
with rough edges is again given by the combination law d2oT ). 
where pj is the scattering matrix of free motion in the region 
ajj-i < x < Xj and Sj is the scattering matrix of the j-th edge 
step. The scattering matrix pj is given by equation d27l i with 
the matrix elements (f28b modified as 



Finally, we follow [17] to specify the scattering matrix Sj. 

Consider the wire shown in figure [5] A single edge step 
at x = divides the wire into the region A (x < 0) and 
region B (x > 0). The widths of the regions A and B are 
Wa = hA — Aa and Wb = hs — ds, respectively. We 
consider the case Wa > Wb and we assume, that the wire 
cross-section Wa includes the wire cross-section Wb (see the 
figure). In this situation, the confining potential is simply 

!0, x < 0, dA < y < hA 
0, x > 0, d B <y<h B (40) 
oo, elsewhere. 

So we can use ([35} and write the wave function at x — as 

N A 

V (0-e,y) = J2 [A+(0)+A-(0)] X *(y), e -> 0, (41) 



Nb 

¥3 (0 + e,y) = ^[B+(Q)+S-(0)]xf(y) I e -J- 0, (42) 

n=l 

where A±(x) = a±e ±ik - x , B±(x) = b±e ±ik * x , and N A 
and Nb are the channel numbers in the regions A and B. The 
continuity equation tp(0 — e, y) = ip(Q + e, y) takes the form 

A+ +A- = C{B -> A){B + +B~), (43) 

where C(B — > A) is the matrix with the matrix elements 



C(B -> A) 



flA 

n.„ ■ i Xm(v)Xn(y) d y 

d A 



(44) 



and dimensions Na X Nb- Similarly, the continuity equation 
§f (0 — e, y) = §f (0 + e, y) can be written in the form 

K A (A + -A~) = C{B -> A)K B {B + -B-) 



(45) 



where K A and K B are matrices with the matrix elements 



4\. 



Ah-* Ax s 



(39) 



(K A ) 



rnn — u mn j 



\ - k B i) 



(46) 
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and dimensions Na X Na and Nb X iVs, respectively. Com- 
bining (l43T > and d45b one finds the matrix equation d2"TT i with 
the scattering matrix Sj = S composed of the matrices 



where 



t = -2MH BA , 

r = -2C(B -» A)MH BA - I A , 
t' = C(B -> A)[MN + I B ], 
r' = MN, 



H BA = -{K B )- l C T (B -> ,4)if A , 
N = I B +H BA C(B ^ A), 



(47) 



(48) 



with / and I B being the unit matrices and C T being the 
matrix obtained by transposition of the matrix C. The dimen- 
sions of the matrices t, r, t' and r' are Nb x A/^, x Na, 
Na x Nb and Nb x A^b, respectively. 

Proceeding in a similar way one can derive sj for the situ- 
ation Wa < Wb, with the cross-section W A included in the 
cross-section W B . In this case Sj is composed of the matrices 



t = C(A -4 S)[MAT + / A ], 

r = MN, 

t' = -2MH A b, 

r' = -2C(A -> B)MH A b - I B , 



(49) 



where Af, Hab-, and N are the matrices (08]) with the index 
yl replaced by B and vice versa. 



C. Averaging over the samples made of the building blocks 

The conductance (JTTJ needs to be evaluated for a large en- 
semble of macroscopically identical wires because it fluctu- 
ates from wire to wire. To evaluate the ensemble-averaged 
results like (g), (g 2 ) — (g) 2 , (l/g), etc., we need to perform 
averaging typically over 10 3 — 10 4 samples with different mi- 
croscopic configurations of disorder. Moreover, the ensemble 
averages are studied in dependence on the wire length. Espe- 
cially for long wires (L ~ £) with a large number of conduct- 
ing channels (N c ~ 30 — 300) already the scattering-matrix 
calculation of a single disordered sample takes a lot of com- 
putational time. To decrease the total computational time sub- 
stantially (say a few orders of magnitude), we introduce a few 
efficient tricks, partly motivated by the work 13211 . 

We recall that two scattering matrices, Si and S2, can be 
combined by means of the operation (1241 . written symboli- 
cally as S\ ® S2. This operation is associative, i.e., 



{Si ® S 2 ) ® S 3 = Si ® (S 2 ® S 3 ). 



(50) 



The property (I50t allows us to proceed as follows. We can 
construct the disordered wire of length L by joining a large 
number of short wires (building blocks), where each block 
has the same length (Lb) and contains the same number (rib) 



of scatterers (impurities or edge steps). The scattering matrix 
of such wire can be expressed as 



5 = h O b 2 ® • • • = (pi ® si ® 

(Pn b + 1 <8> S rlb + i (g) ... S 2 n b ) 



,) 



(51) 



where 61 = pi®sit 



is the scattering matrix of the first 



block, 62 = Pn b +i<S>s nb +i® ■ ■ ■ S2n b is the scattering matrix of 
the second block, etc. In the simulation, we can create the Nb 
different blocks (typically Nb = 100) so that we select at ran- 
dom the positions of the rib scatterers in a given block. Then 
we evaluate the scattering matrices (bi) of all Nb blocks. We 
can now readily study the wire lengths L — Lb, 2Lb, iLb, ■ ■ ■ 
by applying one of the two approaches described below. 

A single sample of length L = jLb can be constructed by 
joining j blocks, where each block is chosen at random from 
the Nf, blocks. Clearly, one can in principle construct (Nb) J 
different samples of length jLb, but in practice a much smaller 
number of samples (~ 10 3 — 10 4 ) is sufficient. The 5-matrix 
of each sample can be evaluated by means of ( l5TT l and en- 
semble averaging can be performed. Already this approach 
works much faster than the approach which combines in each 
sample the S-matrices of all individual scatterers. A further 
significant improvement is achieved as follows. 

As before, we evaluate the scattering matrices bi for all 
Nb blocks, but we apply a more sophisticated algorithm: (i) 
We choose at random a single fri-matrix describing a spe- 
cific sample of length L = Lb- (ii) Choosing at random an- 
other bi and combining it with the previous one by means of 
(|5TT i we obtain the 5-matrix of a specific sample of length 
L = 2Lb- (iii) Choosing at random another bi and combining 
it with the 5 matrix obtained in the preceding step we ob- 
tain the 5-matrix of a specific sample of length L = 3Lb- 
In this way we obtain a set of the Landauer conductances 
{G(jLb)}j = i t 2.... for a set of the specific samples with lengths 
L = Lb, 2Lb, .... Repeating the algorithm again we obtain a 
new set {G(jLb)} 3=1,2,...- Repeating the algorithm say 10 3 
times we obtain 10 3 various sets of {G(jLb)}j = i t 2.... and we 
perform ensemble averaging separately for each j. This ap- 
proach saves a lot of time because the 5 matrix of the sample 
of length L = jLb is created by combining a single 6;-matrix 
with the 5 matrix of a sample of length L = (j — \)Lb- 

In reality the positions of all scatterers differ from sam- 
ple to sample. If we take this into account in our simulation, 
the ensemble-averaged results are the same (within statistical 
noise) as those obtained by means of our trick, but the ensem- 
ble averaging takes far much computational time. Owing to 
our trick we can analyze much larger systems, which is essen- 
tial for a successful observation of our major results. 



III. SEMICLASSICAL CONDUCTIVITY OF THE Q1D WIRE 

In this section, the semiclassical Q1D conductivity expres- 
sions are derived both for the impurity disorder and edge 
roughness. Precisely, we assume that the electron motion is 



7 



semiclassical along the direction parallel with the wire and 
quantized in the perpendicular direction. Our approach is 
technically similar to the previous studies of the Q1D wires 
lEEB and Q2D slabs SIM HQ. 

The semiclassical conductivity of the Q1D wire is given as 



FW 



EE 



hk 



LJ m 



fn(k), 



(52) 



where F is the electric field applied along the wire, f n (k) is 
the electron distribution function in the ?i-th conducting chan- 
nel and the factor of 2 includes two spin orientations. Simi- 
larly to the usual textbook approach, we express f n (k) as 

/„(*) = f(E nk ) + ^ df f k nk) r n (E nk ), (53) 

where f{E nk ) = 1/ (exp[{E nk - E F )/k B T] + 1) is the equi- 
librium occupation number of the electron state with energy 
E n k — e« + h 2 k 2 /2m and T n {E nk ) is the relaxation time. 
Setting (l53l into d52l we obtain at zero temperature 



2e 2 



JV 



irmW 



E k n T n (E F ), 



(54) 



where k n is the Fermi wave vector in the channel n [equation 
[8). The function (l53l obeys the linearized Boltzmann equation 



eF df(E nh 
h dk 



= EE^»'( fc > fc ') 

n' k' 
X [fn>(k')-fn(k)}, 



(55) 



where 

W n<n ,%k') 



2tt 



i',k'\U\n,k)\l v S(E nk -E n , k ,) (56) 



is the Fermi-golden-rule probability of scattering from \n, k) 
to \n', k'), with U being the scattering-perturbation potential 
and \n,k) = L^ 1 / 2 exp(ikx)xn(y) being the unperturbed 
electron state. The index av means that W n , n ' (k, k') is av- 
eraged over different configurations of disorder. 

From (|55l l we find (see the appendix A) the relaxation time 



T n {E F ) = 



m \ - 



2TT 2 h^ 

n' 

where K is the matrix with matrix elements 

6 nn > ^ \(n, k\U\n, k')\ 



(57) 



K nn , = ^EE 

k k' 



r\\2 
av 



xk 2 S{E nk - E F )S{E llkf -E F )-\(n,k\U\n',k')\ 



f\\2 
av 



xkk'S(E nk ~ E F )S(E ri 



E F ) 



(58) 



The Boltzmann Q1D conductivity thus reads 



2e 2 1 



EE*»M^- (59) 



n— 1 n' — l 



A. Semiclassical conductivity of the Q1D wire with impurities 

If we set for U the impurity potential ©, the matrix ele- 
ments d58l can be expressed (see the appendix B) in the form 



K , — 




where 



7 = m^/H 2 



Using d60b we obtain from d59l > the Q1D conductivity 



2e 2 1 ^ kl 



h ^ 2 rii ^ 1 + tp N ° 



fc„ 



(60) 



(61) 



(62) 



For N c 1 the sum in d62l converges Xokp/2 and d62l con- 
verges to the 2D limit 



2e 2 fc F fc F 
h 2 ^ 2 ni ' 

derivable from the 2D Boltzmann equation. 



(63) 



B. Semiclassical conductivity of the Q1D wire with rough edges 



To evaluate the conductivity ( 1591 ) for the wire with rough 
edges, we need to determine the perturbation potential U pro- 
duced by the edge roughness potential in the figure Q] to set 
the resulting U into the right hand side of d58l) . and to evaluate 
K nn /. All this is performed in the appendix C. The result is 



7T 2 ^ 2 

w 6 



No 



E»Vr[^-*»l)- 



n 2 n [F(\k n - k n >\)- 



•F(\k n ~\~ k n i | 



(64) 



where 8 is the root mean square of the fluctuations of the edge 
coordinates d(x) and h(x) (see below) and F(q) is the Fourier 
transform of the roughness correlation function F(x). 

To specify 6, F(x), and F{q) we recall (see figure Q]), 
that d(x) and h(x) fluctuate with varying x in the intervals 
(—A, A) and (W — A, W + A), respectively, by changing 
their values abruptly after constant steps Ax. Obviously, the 
values of d(x) are distributed in the interval (—A, A) with the 
box-shaped distribution. Using such distribution we find, that 



(d(x) 2 ) - (d(x)) 2 = (d(x) 2 ) = A 2 /3 = S 2 



(65) 



where the symbol S labels the root mean square of d(x). The 
correlation function F(x) is defined as 



(d(x! +x)d(x 1 )) = S 2 F(x). 



(66) 
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FIG. 6: The mean resistance (p) as a function of the wire length L. 
Note that (p) is reduced by the contact resistance p c = 1/N C and 
L is scaled by the mean free path /. The left panel shows the mean 
resistance of the wire with impurity disorder, the right one shows the 
mean resistance of the wire with rough edges. The wire width W 
is fixed to the value W/\f = 17.4, which means that N c = 34. 
The parameters of disorder are listed in the figure together with the 
resulting values of I. Determination of I is demonstrated in the next 
figure. Note that for the Au wire (Ef = 5.6eV and m = 9.109 x 
10~ 31 kg) we have \ F = 0.52nm and W = 17.3A F = 9nm, with 
the smallest A being only O.Olnm. Such small A is obviously not 
realistic, but it is used to emulate the weak roughness limit. 



impurities edge roughness 




Lll Lll 



FIG. 7: The full lines show the selected numerical data from the pre- 
ceding figure. The dashed line in the left panel shows the linear fit 
(p) = p c + pdifL/W, while the dashed line in the right panel is 
the linear fit (p) — p%'' + pdifL/W, with p%' ' being the effective 
contact resistance. In both cases the resistivity pdi / is a fitting param- 
eter and I is extracted from the Drude formula pdif = 2/(kpl)- To 
determine p e J^ , the dashed line in the right panel is extrapolated to 
L — 0. We find p e J 1 ~ 1/7.5 while p c = 1 /34. Inset shows the full 
line from the right panel for L/i <g 1, where the transport is quasi- 
ballistic. The dotted lines show the linear fit (p) — p c + p q bL/W, 
where the quasi-ballistic resistivity p q t is a fitting parameter. The 
quasi-ballistic mean free path l q t, is given as l q t = 2/(kpp q b)- 



We use the box distribution and we take into account that the 
step Ax plays the role of the correlation length. We obtain 



F(x) = 



1 - \x\/Ax, \x\ < Ax 



0, \x\ > Ax 

The Fourier transform of the last equation is 

2(1 - cos(Axq)) 



T(q) = Ax- 



(Axq) 2 



(67) 



(68) 



The same results as for d(x) hold also for h(x) — W, because 
the roughness of both edges is the same. 



IV. NUMERICAL AND ANALYTICAL RESULTS 

In subsection A, the quantum transport in the wires with im- 
purity disorder and wires with rough edges is simulated in the 
quasi-ballistic, diffusive, and localized regimes. In subsection 
B the crossover from the quasi-ballistic to diffusive regime is 
explained by means of an intuitive model based on the concept 
of the open channels. A microscopic analytical theory of the 
crossover is given in subsection C. In subsection D, the diffu- 
sive mean-free path obtained from the Landauer conductance 
and mean-free path from the Boltzmann theory are compared 
for both types of disorder. In subsection E the wire with rough 
edges is studied for large roughness-correlation lengths. 

In principle, all our transport results can be expressed 
and presented in dependence on the dimensionless variables 
A/Xp, Ax/Xp, W/Xp, etc., with the Fermi wave length Xp 
being the length unit. Nevertheless, in a few cases we also use 



the normal (not dimensionless) variables and we present the 
results for the material parameters m = 9.109 x 10~ 31 kg and 
Ep = 5.6eV (Xp = 0.52nm), typical of the Au wires. 



A. Quantum transport in wires with impurities and rough edges 

We mostly simulate the wires with the number of the con- 
ducting channels being N c = 34. This number emulates the 
limit iV c > 1 without spending too much computational time. 
Whenever needed we also use much larger N c , the largest one 
being N c — 347. We calculate the Landauer conductance for 
10 4 wires and we evaluate the means. 

We start with discussion of the mean resistance (p). In the 
figure|6]the mean resistance of the wire with impurity disorder 
is compared with the mean resistance of the wire with rough 
edges for various parameters of disorder. The data obtained 
for various parameters tend to collapse to a single curve when 
plotted in dependence on the ratio L/l. Hence it is sufficient 
to discuss only the data for one specific choice of parameters. 

The figure [7] shows the selected numerical data (full lines) 
from the figure [6] In accord with the textbooks OJ-Slj the 
mean resistance of the wire with impurities follows for L > 
the standard diffusive dependence (p) = p c + pdifL/W. 
However, the mean resistance of the wire with rough edges 
shows a more complex behavior. For L — > it follows the 
linear dependence (p) = p c + p q ^L/W, where p q b is the 
quasi-ballistic resistivity. Only for large enough L it shows 
crossover to the diffusive dependence (p) — p e J' +pdifL/W, 
where the resistivity p^if is much smaller that p q \, and the 
effective contact resistance p e J* strongly exceeds the funda- 
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FIG. 8: Typical conductance (hip), mean resistance (p) and mean 
conductance (g) versus !//£. The results for the wire with impurity 
disorder (left panels) are compared with the results for the wire with 
rough edges (right panels). The calculations were performed for var- 
ious sets of the parameters, listed in the figure [6] The results for 
various parameter sets collapse to a single curve shown in a full line 
(see the remarks in the text). The dashed lines in the top panels show 
the fit (ing) = —L/£, from which we determine the localization 
length £. The resulting values of £ are shown. 



mental contact resistance p c . In other words, the wire with 
rough edges shows two different linear regimes (the quasi- 
ballistic one and the diffusive one) separated by crossover, 
while the wire with impurities shows a single linear regime 
for the quasi-ballistic as well as diffusive transport. 

The figure|7]also shows that the mean resistance of the wire 
with impurities increases for L/l > 10 slightly faster than 
linearly, which is due to the weak localization. However, the 
mean resistance of the wire with rough edges increases lin- 
early even for 10 < L/l < 30. The origin of this difference 
will become clear soon. 

The figure [8] shows the numerical results (full lines) for the 
the typical conductance (In g), mean resistance (p), and mean 
conductance (g) in dependence on the ratio L/£. The cal- 
culations were performed for various sets of the parameters, 
shown in the figure [6] The results for various sets tend to 
collapse to a single curve (full line) when plotted in depen- 
dence on L/£. We see for both types of disorder, that the 
numerical data for (In g) approach at large L the dependence 
(lng) = This is a sign of the localization 

Fitting of the numerical data provides the values of £ shown 
in the figure. In the wire with impurities we find the result 
£/l ~ 0.9A^ C , which agrees with the theoretical [37] predic- 
tion ^/l — N c and with numerical studies [6]. In the wire with 
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FIG. 9: The full lines show our numerical data for the mean resis- 
tance (p), inverse mean conductance l/{g), and conductance fluc- 
tuations y/vstr(g). The dashed lines in the top panels are the linear 
fits (p) = p c + p d ifL/W (left panel) and (p) = p% n + pdifL/W 
(right panel). The dashed lines in the middle panels are the linear fits 
V<S> =Pc + PdifL/W (left panel) andl/O?) = p%" + p m L/W 
(right panel). Onset of strong localization is marked by arrows at the 
points, where the numerical data start to deviate from the linear fit 
remarkably. The relative deviation from the linear fit is shown in a 
separate figure [figure dlOH. The dotted line in the bottom panels 
shows the theoretical value (0.365) of the conductance fluctuations, 
predicted in the limit £/£ <C L/£ <C 1 for the white-noise disorder. 



rough edges we find £/Z ~ 1.4iV c . This does not contradict 
26ll . which reports £/1±d — N c , but lip is ir/2 



the work [21 



times larger than our I. Finally, due to the localization also (p) 
and (g) depend on L/£ exponentially at large L/£. 

The figure|9]show the mean resistance (p) and inverse mean 
conductance 1 /(g), now in the linear scale. Concerning the 
mean resistance, the wire with rough edges and wire with im- 
purities behave similarly: (p) rises with L/£ linearly on the 
scale Z/£ <C L/£ < 1, as is typical for the diffusive regime. 
However, both types of wires show a quite different 1 / (g). In 
the wire with impurities 1/ (g) rises with L/£ linearly in the 
interval Z/£ -C L/£ < 1, while in the wire with rough edges 
1/ (g) shows the linear rise with L/£ in the interval as large 
as Z/£ -C L/£ < 2. The (g) oc 1/L dependence is a sign 
of the diffusive conductance regime, which now persists up to 
L/£ ~ 2 and which was not observed in lf20h due to the too 
narrow length window (as explained by the authors). 

Further, the slope of the dashed lines is larger for the edge 
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FIG. 10: Numerical data for (p) and 1 / (g) from the preceding figure, 
presented as a relative deviation from the linear fit. The dashed line 
shows the weak-localization-mediated relative deviation \pdif^, 
shifted by replacement L —¥ (L — 8/) to obey the limit L>l. 



roughness than for the impurities. This can be understood if 
we write the equation (p) ~ l/(g) ~ (2/tt)(€/N c I)(L/£) and 
we realize that the ratio £/N c l is larger for the edge roughness. 

The figure [9] also shows the conductance fluctuations. The 
fluctuations in the wire with impurities approach the universal 
value 0.365, derived HQ in the limit Z/£ < L/£ <C 1 for 
the white-noise disorder. It is remarkable that the fluctuations 
in the wire with rough edges show a length-independent uni- 
versal value (of size ~ 0.3) just in the interval l/^ <C L/£ < 
2, in which we see the linear rise of 1/ (g). Coexistence of 
the universal conductance fluctuations with the conductance 
~ £/L is typical of the diffusive conductance regime J2]. 

In the figure [9] onset of strong localization is visible on the 
first glance at the points (marked by arrows), where the nu- 
merical data for (p) and l/(g) start to deviate from the linear 
dependence remarkably. For the edge roughness the inverse 
conductance shows onset of localization at L ~ 2£ : note that 
the corresponding conductance fluctuations are not universal 
just for L > 2£ (they decay with L). 

The figure [10] shows the relative deviation from the linear 
fit, obtained from the numerical data in figure [9] Also shown 
is the relative deviation \pdif^, obtained from the formula 
( TToT i. As expected, the inverse conductance of the wire with 
impurities exhibits for I <C L < £ the deviation close to 
\pdif~W' ^his * s evidently not the case for the wire with 
rough edges. First, if L < 0.2£, both (p) and I /(g) show 
a large negative deviation due to the crossover from the quasi- 
ballistic to diffusive regime. Second, if 0.2£ < L < 2£, then 
1/ (g) shows the deviation as small as < 0.08 and almost no 
deviation for 0.4£ < L < 1.1£. In other words, l/(<?) ex- 
hibits up to L ~ 2£ the linear diffusive behavior with a minor 
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FIG. 11: Mean resistance (p), inverse mean conductance l/(g), and 
conductance fluctuations ^var(g) for various wire parameters. 



nonlinear deviation. Finally, (p) shows a steeply increasing 
deviation at L ~ £ due to the localization. 

We have sofar discussed the numerical data for the wire pa- 
rameters listed in figure[6] Apart from small differences due to 
the statistical noise, these data collapse almost precisely to the 
same curve, when plotted in dependence on In fact, such 
single-parameter scaling (dependence solely on L/Q holds 
exactly for weak disorder [ 19, 38, 39]. If disorder is not weak, 
the data can deviate from the single-parameter scaling and the 
question is whether our findings hold generally. 

For instance, already in the figure [6] we do not see for var- 
ious parameters exactly the same curves. However, the dif- 
ference between various curves is so small, that the findings 
extracted from one of these curves (see figure [7]) are obtain- 
able with a minor quantitative change from other curves. So 
we expect that the same findings hold for any (reasonable) 
choice of the wire parameters. A strong support for this ex- 
pectation is that the parameters used in the figure [6] are very 
different: the mean free paths I range from I/Xf = 41 up to 
l/X F = 3.6 x 10 6 . 

The figurefTTlshows again the numerical data for (p), 1/ (g), 
and y/var(g), but for various sets of the wire parameters. Ob- 
viously, the data for various sets do not collapse exactly to the 
same curve. This may be due to the fact that disorder is not 
weak, however, the resulting curves are also sensitive to how 
accurately we determine £. We could improve proximity of 
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FIG. 12: Effective number of open channels, N% , versus the 
roughness correlation length Ax for various N c and various rough- 
ness amplitudes A. The parameters of the edge roughness are 
scaled as Ax/\ F and A/W. The results for A/W = 1/18 are 
shown by open circles, the dashed and full lines show the results for 
A/W = 1/180 and A/W = 1/900, respectively. 



the curves in the figure [TT]by simulating a larger ensemble of 
samples and a larger wire length (in order to obtain a more 
accurate £). However, the presented proximity is quite suffi- 
cient in the sense that each of the curves allows to obtain the 
results very similar to those in figures [9] and [10] Proximity of 
the curves is satisfactory also with regards to the fact that the 
values of £ obtained for various parameters in the figure [TT] 
vary in the range of five orders of magnitude. In this respect 
we can also say, that the figure[TT]confirms universality of the 
conductance fluctuations in the wire with rough edges: they 
are of size y/var(g) ~ 0.3, reported by others II 9l l25l 1271 12811 . 



B. Crossover from the quasi-ballistic to diffusive transport in 
wires with rough edges: Intuitive analytical derivation 

Let us examine the crossover from (p) — p c + p q \,L/W 
to (p) = p e Jf + pdifL/W, observed in the figure 111 Such 



crossover was not observed in the works 11211 126H . where a 
similar situation was studied numerically. Therefore, we first 
analyze the conditions of observability. We define the effec- 
tive number of the open channels, Af?" = 1/ p%'* ', and we 
evaluate N%'* numerically (by means of the same procedure 
as in the figure |7]l for various wire parameters. 

The figure[T2]shows JV® in dependence on the roughness- 
correlation length Ax for various N c and various A. We see 
that N%** reaches for Ax/\p — > a minimum value which 
is roughly 6 and which depends, within our numerical accu- 
racy, neither on N c nor on A/W. In other words, N*" ap- 
proaches for Ax — > the universal value ~ 6. The calcula- 
tions in i2ll l26ll were performed for Ax ~ O.BXp, but only 
for N c = 5. This value is obviously too small for noticing the 



impurities 



edge roughness 




A 
V 



10 



L/£ = 0.5 

^0000000000000000000000000000000 



10"' 



10 



e>/v L/§ = 0.5 



'"••oooooo, 







10 20 30 
n 



10 20 30 
n 



FIG. 13: The top panels show the transmission probability (T n ) ver- 
sus L/£ for the channel indices n = 1,2, ... N c , where N c = 34. 
For n ordered increasingly, the resulting curves are ordered decreas- 
ingly: the top curve shows (T„-i), the bottom one shows (T n =N c ). 
The bottom panels show (T n ) versus n for various L/£. 



existence of Nj?ff ~ 6. Hence the work [21, 26] reported the 
diffusive dependence (1/g) ~ 1/N C + pdifL/W rather than 
the dependence (1/g) ~ 1/A r c e// + p dl fL/W. Finally, for 
Ax/\f ^> 1 we see that N%f* approaches N c . That limit is 
studied in the last subsection. 

Further, we look at the numerical data for (T n ), shown in 
the figure Q~3] The theory based on the white-noise disorder 
predicts, that the conducting channels are equivalent iTTsL \l$\ 
in the sense that (Ti) = (T 2 ) ■ ■ ■ = (T Nc ). In the figure [13] 
this equivalency is reasonably confirmed for the wire with im- 
purity disorder but not for the wire with rough edges. In the 
wire with rough edges (T n ) decays fast with rasing n which 
is easy to understand classically: in the channel n = 1 the 
electron avoids the edges by moving in parallel with them, 
while in the channel n = N c the motion is almost perpendic- 
ular to the edges, resulting in frequent collisions with them. 
As a result, the (T n ) dependence in the right panels of fig- 
ure [13] shows for L ~ £ the coexistence of the quasi-ballistic, 
diffusive, and strongly-localized channels, already reported in 
previous works JH |2(1 . 

Concerning the coexistence, two comments are needed. Ev- 
idently, the coexistence is not in contradiction with the fact 
that all T n decay in semilogarthmic scale linearly with a sin- 
gle parameter L/f, when L > ( (see also [8]). Further, it 
is clear that if the diffusive regime means T n oc 1/L in all 
N c channels at the same value of L, then there is no diffusive 
regime but only crossover from the quasi-ballistic regime to 
the localization regime J20ll . In the preceding text we were 
speaking about the diffusive regime in the sense J2II l26Tl that 
(p) oc L and I /(g) oc L. 
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For the uncorrelated impurity disorder, the transmission 
(T n ) in absence of the wave interference takes the form J^] 



(T n ) =L n /(L n + L), 



(69) 



where L n is the characteristic length. For the edge roughness 
we can adopt ( |69l as an ansatz. We will prove later on that 
the ansatz is indeed correct for the uncorrelated roughness. In 
what follows we combine the ansatz (169b with the concept of 
the open channels and we explain all major features of the 
crossover from the quasi-ballistic regime to the diffusive one, 
albeit the formula d69t cannot capture the fact that besides 
the quasi-ballistic and diffusive channels there are also the 
localized ones. 

We asses L n from the numerical data in figure [13] From 
the figure[[3]it is obvious, that L n in the channels with n ^> 1 
is much smaller than L n in the channels with n — ► 1. We 
emulate these findings by a simple model. We introduce 
the number <C N c and we assume, that the channels 

n = 1,2,..., iVg*' have the characteristic length L a , while 
the rest of the channels has the characteristic length Lb <C L a . 

By means of the above model, we can estimate the mean 
resistance (p) and mean conductance (g) for the wire lengths 
< L < £. From the figure|9]we see that 



(p)^l/(g), 0<L<1 
We therefore rely on the equation 
1 



(9) 



En=l (-^") 



■■It 



Y, n l N eff + 1 ( T n) 



(70) 



(71) 



In the quasi-ballistic limit (L <C Lb) the first term in the de- 
nominator of ( TTTb is simply '* and the second term can be 
evaluated by means of d69l ). For L n = Lb we find 



1 N c - N c ^f T 1 
{p)C "T c + ^!Lb- L ^N- c 



1 



NrLh 



(72) 



where the right hand side holds for N^'f <C N c . In the dif- 
fusive regime {I -C L < £) we evaluate the denominator of 
( TTTb by means of d69l and we neglect the second term in the 
denominator assuming that N%'* L a ^> N c Lb- We get 



(P) 



Lg + L 

N e J f L a 



N, 



eff 



N! ff L a 



L, 



(73) 



where the right hand side is the limit L » L a . This means 
that we assume (T n ) ~ L n / Lfor all N^** channels, i.e., we 
ignore that the channel n = 1 is almost quasi-ballistic even at 
L ~ £ (see fisure \13h Nevertheless, we succeed to obtain all 
major features of the crossover from the quasi-ballistic regime 
to the diffusive one (see mainly the next subsection). 

If we compare the formulae d72l > and d73l with the formulae 
(p) = p c + PqbL/W and (p) = p e Jt + p dif L/W, we obtain 



We have assumed above that N^^L a 3> N c Lb- This means 
that p q b 3> Pdif- Indeed, we will see that p q b/ Pdif — N c 



The formula d73l holds only for L < as the equation ( TTOb 
holds for L < However, the formula 

1 1 1 



N, 



If 



m n L 



-L 



(75) 



holds for I -C L < 2^ as we know from our numerical 
data. (The fact that the conductance behaves diffusively up 
to L ~ 2£ has previously been recognized from the conduc- 
tance distribution B22I1 . We do not show here the conductance 
distributions as they are similar to those in B22ll24ll .) 

The formulae d72l i and ( T73b show, that the crossover from 
(ft) = p c + PqbL/W to (p) = p e J f + p dlf L/W is due to 
the channel non-equivalency. In what follows, the formulae 
(|72| > and (l73l will be derived from the first principles, without 
using the parameters L a and Lb- 



C. Crossover from the quasi-ballistic to diffusive transport: 
Microscopic analytical derivation 



We express the transmission probability (T n ) as 



1 



rn — 1 



(76) 



where R mn is the probability that an electron impinging the 
disordered region in the m-th channel is reflected back into the 
ri-th channel. In the work [20] the wire with rough edges was 
analyzed in the quasi-ballistic limit and the reflection proba- 
bility {R mn ) was derived by means of the first order perturba- 
tion theory. The result is 



K m hu n 



T(\k n 



)L. 



(77) 



where K n — (ir/W)n and the factor of 2 accounts for two 
edges. (In fact, the result given in |20ll involves a missprint. 
The result ( |77] i can also be extracted from the backscattering 
length reported in ll40[|4lll .) 

The mean conductance in the quasi-ballistic limit reads 



(9) 
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where l q b is the quasi-ballistic mean free path: 
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(78) 



(79) 



where S 2 — A 2 /3. Note that the quasi-ballistic limit 
contains only the backscattering contribution oc F( \ k n + k m \ ) , 
while in the diffusive regime (equation |64|| also the forward- 
scattering contribution oc J-(\k n — k m \) is present. The mean 
resistance in the quasi-ballistic limit is 



13 



^ 10 

£ 
3 




A = 0.002 nm 



o o o o o o 



OOOOOOO 0000000 




15 20 25 30 
W [nm] 




£ F [eV] 



FIG. 14: Top panel: Quasi-ballistic mean free path l q b in the Au wire 
as a function of the wire width W for various roughness amplitudes 
A. The roughness correlation length is fixed to Ax = 0.125nm, the 
Au material parameters are m = 9.109 x 10~ 31 kg and Ef = 5.6eV. 
The circles show our numerical data and the squares (connected by 
a full line) are the data points obtained from the expression d79t . 
Bottom panel: The same calculations as in the top panel, but W is 
fixed to 9nm and the Fermi energy is varied. To explore the small A 
limit reliably, we have to use the values of A which are too small to 
be realistic. 



In the figure [14] the expression d79l is compared with l q i, 
determined numerically by means of the approach discussed 
in figure [7] (the right panel and inset to the right panel). The 
formula d79j > agrees with our numerical data if the roughness 
amplitude A is small. This is what one expects, because the 
perturbation expression ( F79b is exact in the limit A — > and 
our scattering matrix calculation is (in principle) exact for any 
A. As A increases, the result d79l ) fails to agree with our 
numerical data because the scattering is not weak [42]. 

Simple formulae can be derived for l q i, and I if the rough- 
ness is uncorrected (Ax/Xp <C 1). We start with I. If 
Ax/Xp <C 1, the correlation function ((68) is simply J-(q) = 
Ax. Consequently, the backscattering and forward-scattering 
terms in d64l become the same and the matrix d64l > reduces to 
the diagonal form 
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3W 6 



X £ nV 



(81) 



where the factor of 2 is just due the equal contribution of the 
backward and forward scattering. We set (l8lT l into the Boltz- 
mann conductivity (l59l l and we extract from (l59l the diffusive 
mean free path. It reads 
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(82) 



For N c ^> 1 the summations in 

N c 



— ~ — k F , 
n z 6 

n— 1 
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can be approximated as 

W 3 



p 

4^ F ' 



and the semiclassical diffusive mean-free path becomes 
l 2W 2 



A 2 k 2 F Ax' 



(83) 



(84) 



Now we evaluate for Ax/\p <C 1 the quasi-ballistic mean 
free path. We rewrite d79l i into the form 
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where 

Setting into 

l n 
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the formula T(q) 
3W 2 k m k n 



Ax we obtain 



3W 6 k m k„ 



2A 2 K 2 n K 2 l Ax 2ir 4 A 2 m 2 n 2 Ax' 
Combining (|85l l with (|87] i and using d83l we obtain 
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(85) 
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(88) 



We recall that l q t, is limited exclusively by backscattering. 

The formulae ([84]i and ([88]l hold for Ax/X F <C 1. In figure 
[TBIwe compare them with the original formulae valid for any 
Ax. We can see that the major difference is absence of the 
oscillating behavior in the formulae d84i l and d88b . 

Finally, the ratio l/l q b can be expressed as 



J_ 

lqb 



24 



(89) 



The result ( 1891 is universal - independent on the wire param- 
eters and parameters of disorder. The figure [16] shows that 
the universal relation l/l g b oc N c is confirmed by our exact 
quantum-transport calculation. Obviously, since the formula 
(|84| > is the semiclassical Boltzmann-equation result and for- 
mula d88l is the weak-scattering limit, the formula d89b cannot 
reproduce the exact quantum results quantitatively. 

The diffusive mean-free path ( f84b contains both the back- 
ward and forward scattering. It is therefore interesting, that 
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FIG. 15: Quasi-ballistic mean free path l q t and diffusive mean free 
path / in the Au wire with rough edges, calculated as a function of the 
wire width W and Fermi energy Ef- The thick lines show the for- 
mulae {88} and J84t . derived for the uncorrelated roughness. The thin 
lines show the formulae valid for an arbitrary correlation length Ax, 
namely the quasi-ballistic result ( |79b and the semiclassical / extracted 
from the Boltzmann conductivity l |59t . The roughness amplitude A 
and roughness correlation length Ax are fixed to A = 0.5nm and 
Ax = 0.125nm (the limit Ax/\ F < 1 is fulfilled for all consid- 
ered data). The Fermi energy used in the top panels is Ef = 5.6eV, 
the wire width in the bottom panels is W = 9nm. The oscillations 
with sharp minima appear whenever the Fermi energy approaches 
the bottom of the energy subband n — N c . 



the same result can be obtained when the quasi-ballistic 
(backward-scattering-limited) resistance (fSOb is extrapolated 
into the diffusive regime by means of the ansatz d69l l. We start 
from 

-, -l 



(p) * 1/(9) = 



E 

.n=l 



(T n ) 



(90) 



and we use the ansatz d69j l. In the quasi-ballistic limit (L <C 
L n ) we obtain from d69l the formula (T n ) = 1 — L/L n . We 
set this formula into d90l ) and we compare d90b with the quasi- 
ballistic expressions ( f80b and d85l l. We find that 
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(91) 



In the diffusive limit (L n <C L) we obtain from d69l the for- 
mula (T n ) = L n /L and from d90l > the diffusive expression 
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We compare this expression with p = Pdif-^, where p c uf 
2/(kpl). We find the mean free path 
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FIG. 16: Ratio l/l q b, where / is the diffusive mean free path and l q t 
is the quasi-ballistic mean free path, evaluated as a function of N c for 
various roughness amplitudes A. The roughness correlation length 
Ax is fixed to the value Ax/Xf = 0.24 which reasonably emulates 
the uncorrelated roughness. The open symbols show the results of 
our quantum transport simulation, where / and l q t are calculated as 
in the figure|7] 



If we set into d93l ) the uncorrelated limit d87| i, we obtain again 
the Boltzmann mean-free path d84b . This is the proof that the 
ansatz d69l works correctly for the uncorrelated roughness. 
Now it is useful to make two remarks. 

First, our characteristic length L n should not be confused 
with the often used J4lll43ll attenuation length L n . Our L n is 
defined by the ansatz do*9l and we have just seen, that the ex- 
pression d93l l gives for such L n the mean free path coinciding 
with the mean free path obtained from the Boltzmann equa- 
tion. This is the momentum-relaxation-time-limited mean- 
free path. However, if one sets into d93l l the attenuation length 
(equation (5.2) in l43ll ). one obtains from (f93T > the scattering- 
time-limited mean-free path, i.e., the mean distance between 
two subsequent collisions. For the uncorrelated roughness the 
latter is exactly twice shorter than the former one. 

Second, we note that the ansatz d69l ) does not work for arbi- 
trary roughness. Indeed, we can set into the formula d93l ) the 
more general expression ( f86b and we can compare d93l > with 
the Boltzmann mean free path valid for an arbitrary correla- 
tion length Ax. In such case the formula ((93) fails to repro- 
duce the Boltzmann-equation result. 

Assuming the uncorrelated roughness, we are ready to de- 
rive analytically the effective number of the open channels, 
iVg** . The expression d90t can be formally written as 
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where the symbol Nf!" is defined as 
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(95) 



in order to obtain ( |90b again. We now express the transmission 
(T n ) by means of the formula ( f69b and the mean free path I 
by means of the formula d93t . We obtain 



E 



L n + L 



(96) 



In the diffusive limit (L 3> L n ) we obtain after some algebraic 
manipulations the equation 



N eff 



(97) 



The expression d97) no longer depends on L and evidently 
represents the effective number of the open channels, if the 
resistance ( 194-b is considered in the diffusive limit. We set into 
([97) the formulae d£B and ((87). We obtain 
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For JV C > 1 the first sum in the denominator of 



(98) 



becomes 
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(99) 



and other sums in (l98l are already known [see equations (|83T >1 . 
We arrive to the result 



N eff 
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(100) 



which implies that N%** is a universal number depending nei- 
ther on the roughness amplitude A nor on the number of the 
conducting channels, N c . All this agrees with our quantum- 
transport calculation in the figure Q~2j except that the result 
( 1 1001 ) underestimates the numerical value N*ff — 6 about 
twice. We however recall that the result ( 1 1001 relies on the 
formulae d9D and d87) which are not exact. 



D. Diffusive mean free path: quantum-transport results versus 
the semiclassical Boltzmann results 



In this subsection, the mean free path obtained from the 
semiclassical Boltzmann equation is compared with the mean 
free path obtained from the quantum resistivity gdif- For the 
impurity disorder we see a standard result JaLZI]: the semi- 
classical and quantum mean-free paths coincide for weak im- 
purities. On the contrary, for the edge roughness we find, that 
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FIG. 17: Top panel: Diffusive mean free path / in the Au wire with 
impurity disorder as a function of the wire width W, calculated for 
three different sets of the impurity parameters (listed in the bottom 
panel). The semiclassical mean free path obtained from the Boltz- 
mann Q1D conductivity i62\ is shown in a full line: all three sets 
of parameters are intentionally chosen to give the same semiclassical 
result. The data shown by symbols are the quantum-transport results 
obtained from the quantum resistivity pdif (c.f. figure |7J- Bottom 
panel: The same calculations as in the top panel, but W is fixed to 
9nm and the Fermi energy is varied. 



the quantum mean-free path differs from the semiclassical one 
even if the roughness amplitude A is vanishingly small. 

In the figure [17] the quantum and semiclassical mean free 
paths are compared for the impurity disorder. The semiclas- 
sical data are shown in a full line: three different sets of pa- 
rameters are intentionally chosen to provide the same semi- 
classical result. The full lines exhibit oscillations with sharp 
minima, appearing whenever the Fermi energy approaches the 
bottom of the energy subband n = N c . Evidently, the oscil- 
lating curve intersects the quantum results (the data shown by 
symbols), albeit they are slightly affected by statistical noise. 
Both the semiclassical and quantum results follow the trend 
predicted by the semiclassical 2D limit (equation[63l, namely 

1/2 

that I is proportional to Ep and independent on W. In sum- 
mary, in the Q1D wire with impurity disorder the semiclassi- 
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FIG. 18: The right panel shows the diffusive mean free path I in the 
Au wire with rough edges as a function of the wire width W, cal- 
culated for various roughness amplitudes A. The left panel shows 
the same calculation, but W is fixed to 9nm and the Fermi en- 
ergy is varied. The roughness correlation length is kept at the value 
Ax = 0.125nm, which means that Ax/\f < 1 for most of the pre- 
sented data. The full lines show the semiclassical mean free path ob- 
tained from the Boltzmann Q1D conductivity (equations|59|and|64l. 
The open circles show the quantum mean free path extracted from the 
quantum resistivity pdif (c.f. figure [7}. The open squares show the 
mean free path obtained by means of the classical scattering-matrix 
calculation (see the text), in which the localization is absent. To ex- 
plore the small A limit reliably, we have to use the values of A which 
are too small to be realistic. 



cal and quantum mean-free paths coincide. This is a standard 
result [6, 7], known from the theory based on the white-noise 
disorder y|] (see however 14411 ). 

As shown in the figure [18] in the wire with rough edges 
the situation is different. We again compare the semiclassi- 
cal mean free path (full lines) with the quantum mean free 
path (open circles). As before, the full lines exhibit oscilla- 
tions with sharp minima, appearing whenever the Fermi en- 
ergy approaches the bottom of the energy subband n = N c . 
Evidently, the full line and open circles show a quite different 
behavior if the roughness amplitude A is large (notice the re- 
sults for A = 0.5nm). It is tempting to ascribe this difference 
to the weak-perturbation approximation involved in the Boltz- 
mann equation, and similarly, it is tempting to expect that the 
full line will intersect the open circles for sufficiently small A. 
However, this is not the case: even for the smallest considered 
A the open circles are systematically a factor of ~ 2 below the 
full line. In summary, in the wire with rough edges the quan- 
tum mean-free path differs from the semiclassical one (by the 
factor of ~ 2) even if A is vanishingly small. 



To understand the origin of this difference, we now exclude 
from our quantum-transport calculation the wave interference. 
We recall (see subsection II. C), that in the quantum-transport 
calculation the total S'-matrix of the disordered wire is ob- 
tained by combining at random the scattering matrices of 
the building blocks, where the 2N x 2N matrix 6j is com- 
posed of the complex amplitudes t mn ,t' mn ,r mn , and r'. 
To exclude the wave interference, we proceed as follows [7]. 
First, we consider only the conducting channels and we com- 
pletely neglect the evanescent ones: this reduces the size of 
the matrix b$ to 2N C x 2N C . Second, instead of the complex 
matrix 6, we use the real one, in which the complex ampli- 
tudes t ran ,t rnn: v rnn , and 

r 'mn are replaced by the real prob- 



abilities T„ 

R! = \r' 



T' 



and 



\r' mn , respectively. Of course, the wave interference 
is excluded completely, if the length of the building block, if,, 
coincides with the length of the edge step, Ax. Fortunately, 
in practical calculations the wave interference is negligible al- 
ready for Li, ~ I. If we combine the resulting real matrices 
hi by means of the same combination law as before (equation 
l24l >. we obtain the classical transmission probability %nn and 
eventually the classical Landauer conductance 
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Finally, we perform ensemble averaging over many samples. 

In the figure [19] the classical scattering-matrix calculation 
is compared with the quantum one. The mean resistance due 
to the quantum calculation is labeled as (p)n M to distinguish 
from the classical (p) clas - It can be seen that both (p)n M 
and (p) clas exhibit the crossover to the linear diffusive depen- 
dence (p) — p e J* + pdifL/W, but the classical result shows 
a smaller value of pdif and a larger value of p e J'. The value 
p e /f = ~ 1/3.5 is already close to our theoretical 

result N^ff ~ 2.5 (equation 1 1001 ). where the wave interfer- 
ence is excluded as well. The smaller value of p^ t is due to 
the absence of the localization. 

Let us return to the figure [18] The squares show the mean 
free path extracted from (p) c ; as - These data overestimate the 
quantum results (circles) systematically by a factor of ~ 2. 
Compare now the classical scattering-matrix results (squares) 
with the semiclassical Boltzmann results (full lines). The full 
lines intersect the squares only if A is vanishingly small. This 
is due to the fact that the Boltzmann results rely on the weak- 
perturbation theory while the classical scattering-matrix cal- 
culation still involves the exact quantum scattering by a single 
edge step. 

The roughness-limited conductivity derived from the Boltz- 
mann equation (usually further simplified by using the Fuchs- 
Sondheimer model [45, 46]) is often compared with the trans 



port measurements of the metallic nanowires How- 
ever, such comparison cannot verify the Boltzmann result 
unambiguously due to the presence of other processes (like 
the scattering by impurities and grain boundaries) and also 
due to the fact, that the roughness parameters A and Ax 
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FIG. 19: Mean resistance (p) of the wire with rough edges in depen- 
dence on the ratio L/£. The thin full line is the result of the quantum 
scattering-matrix calculation (the same data as in the top right panel 
of figure^, now labeled as (p)q M - The thick full line is the result 
of the classical scattering-matrix calculation, labeled as (n) , . The 

© ' \"f clas 

dashed lines show the linear fits (p) = p e J 1 + pdifL /W , from which 
we determine the effective contact resistance p'g* and the diffusive 
mean free path I = 2/(kFpdif)- The resulting p%'* are shown in 
the figure. Inset shows the classical result again, but for larger L/£ 
in order to stress the linear raise and absence of the localization. 



are not known apriori. Our quantum calculation shows un- 
ambiguously, that the semiclassical Boltzmann approach de- 
scribes the wires with rough edges reasonably, only if A is 
too small to be realistic (in the Au nanowires). In practice, the 
metallic nano- wires are usually fabricated by advanced lift-off 
techniques lf50l - l52tl . which hardly allow to suppress the edge 
roughness to the value A ~ lnm. 

In the table [I] the quantum and semiclassical mean-free 
paths are compared for a more realistic A and W like in the 
figure [18] We see that the quantum result can exceed the 
semiclassical Boltzmann result one order of magnitude. One 
might argue that the quantum result holds only in the coherent 
regime while decoherence is present at any nonzero temper- 
ature. We expect that decoherence will tend to destroy the 
wave interference and to restore the resistance (p) clas result- 
ing from the classical scattering-matrix approach. The result- 
ing classical mean-free path would be twice larger than the 
quantum one [see figureQjO, i.e., it would exceeds the Boltz- 
mann mean-free path in the table|I]by another factor of two. 



E. Wires with strongly correlated edge roughness: Ax/\f 1 

The wires with rough edges have sofar been discussed for 
the roughness-correlation lengths Ax < \p. The only ex- 
ception was the figure [12] where the effective number of the 
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TABLE I: The mean-free path I in the Au wire obtained from the 
quantum-transport simulation is compared with the semiclassical 
Boltzmann mean-free path IboUz (equations 1591 and I64t for various 
values of the wire width W, roughness amplitude A and roughness- 
correlation length Ax. Each set of the parameters considered in 
the table represents strong disorder, for which Isoitz strongly differs 
from I. 



open channels, N%*' , was calculated from Ax C Ap up 
to Ax ^> Xp. We have seen, that the value of N%** ap- 
proaches for Ax ;g> Xp the value of N c , which suggests 
that the diffusive dependence (p) = ljN^** + pdifL/W ap- 
proaches the standard form (p) — 1/N C + pdifL/W, seen in 
the wire with impurities. In this subsection we focus on the 
limit Ax SS> X F . 

The figure [20] shows the same calculations as the figure [9] 
except that the wire with rough edges is now simulated in the 
limit Ax ;g> X F . Unlike the figure [9] the results for the wire 
with rough edges are now close to the results for the wire with 
impurity disorder. 

First, the localization length in the wire with rough edges 
fulfills the relation £/Z ~ 0.92Af c , which is close to the re- 
lation £/l ~ 0.88iV c found for the impurity disorder. This 
is why the resistance (p) and inverse conductance 1 / (g) in 
the figure [20] show for both types of disorder the same lin- 
ear slope with L/£. Moreover, the same linear slope means a 
single linear regime for the quasi-ballistic as well as diffusive 
transport, i.e., the crossover between the quasi-ballistic and 
diffusive regimes, observed for Ax < Xf, tends to disappear. 
(In the figure [9] a remarkably larger slope is seen for the edge 
roughness because £/Z ~ 1AN C .) 

Second, the inverse conductance in the figure [20] shows for 
both types of disorder the linear diffusive regime for L/£ < 1 
and onset of localization for L/£ ~ 1. (In the figure [9] the 
inverse conductance of the wire with rough edges shows onset 
of localization for L/£ ~ 2.) 

Third, the wire with rough edges now exhibits essentially 
the same universal conductance fluctuations as the wire with 
impurities. The size of these fluctuations is in accord with the 
value yj var(g) = 0.365, predicted for the white-noise-like 
disorder. To see such conductance fluctuations for the edge 
roughness is surprising: sofar only the value y/\ds(g) ~ 0.3 
has been reported Ir9 , 25 , 27 , 28 1 . 



To provide insight, it is useful to look at the transmission 
probabilities (T n ). The figure I2TI shows the same calculation 
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FIG. 20: The same calculation and the same symbols as in the 
figure [9] except that the edge roughness is considered in the limit 
Ax/Xf 3§> 1. Specifically, the data obtained for the wire with im- 
purities (left column) are the same as those in the left column of 
figure [9] while for the edge roughness (right column) we show the 
data obtained for Ax/Xf = 19300 and A/Xf = 0.096 (the data 
for A/Xf = 0.019 and A/Xf = 0.96, not shown, look similarly 
except for statistical noise). We recall that the localization length £ 
is determined from the dependence (In g) = —L/£. For the edge 
roughness we now obtain £/Z ~ 0.92iV c , which is close to the impu- 
rity case: f/Z ~ 0.88iV c 



as the figureQj] but with the edge roughness considered in the 
limit Ax/Xf 3^ 1. It can be seen that the figure I2T1 differs 
from the figure Qj] in the following respect: except for a few 
channels with lowest values of n, the wire with rough edges 
and wire with impurities show for the rest of the channels a 
very similar (T„), which is clearly not the case in figure [T3l 
This similarity is responsible for similarity of the transport re- 
sults in the figure[20] The observed similarity could be further 
improved by simulating the wires with N c larger than 34, but 
this is beyond the scope of this work. 

Finally, the figure 1221 shows for Ax ;§> Xp, what happens 
in the wire with rough edges with the diffusive mean free path. 
Specifically, the semiclassical Boltzmann result (full lines) is 
compared with the quantum (open circles) and classical (open 
squares) scattering-matrix calculations. In the limit Ax ^> 
A f all three calculations tend to give the same mean free path. 
This means that the interference effects observed for Ax < 
Xf are no longer important, or in other words, the diffusive 
transport in the wire with rough edges is semiclassical. This 
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FIG. 21: The same calculation and the same symbols as in the fig- 
ure [T3] except that the edge roughness is considered in the limit 
Ax/Xf 3> 1 [we show the data for Ax/Xf = 19300]. 



is again similar to the wire with impurities. 

For Ax ;§> Xp it is easy to show analytically, that the 
semiclassical Boltzmann mean free path rises with Ax lin- 
early. The correlation function F(q) contains the function 
cos(Axg). This function oscillates with q at random and the 
oscillations become very fast for Ax ^> Xp. Consequently, 
the correlation function J-(q) oscillates fast around the value 
F(q) = «L If we ignore these fast oscillations, the diffu- 



sive mean free path can be written as 



I — ~ a a o ^ ' ^ ' k n k n i 



/v F 7T 4 A 2 



(/C) 



1 



(102) 



where 



El 2 K n 
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n 



- (1 - 5 nn >)n 2 n' 2 



1 1 



(k n i fcn) 2 (k n * + fcn) 2 



(103) 



It can be seen that I oc Ax. The I oc Ax dependence resem- 
bles the wire with impurities, where I rises linearly with the 
average distance between the impurities. The formula ( 1102b is 
plotted in the figure [22] in a dashed line. It indeed agrees with 
the numerical data for Ax ^> Xp- On the other hand, the 
formula d84l i. derived in the opposite limit, reasonably agrees 
with the numerical data for Aa; -C Xf, but the dependence 
I oc 1/Ax has no similarity with the impurity case. 

Finally, as also pointed out in the figure caption, the figure 
1221 shows that the semiclassical Boltzmann result (full line) re- 
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FIG. 22: Diffusive mean-free path / in the wire with rough edges as 
a function of the roughness-correlation length Aa; for various rough- 
ness amplitudes A, with I, Ax and A scaled by Af- The wire 
width W is fixed to the value W/Af = 17.4, which means that 
N c = 34). The full lines show the semiclassical mean-free path ob- 
tained from the Boltzmann conductivity (equations |59|and[64l. The 
dashed lines for Ax/Af < 1 show the uncorrelated limit d84t , the 
dashed lines for Ax/Af > 1 show the correlated limit <102t . The 
circles show the quantum mean-free path extracted from the quan- 
tum resistivity pdif (c.f. figure |7](. The squares show the mean-free 
path obtained by means of the classical scattering-matrix calculation 
in which the localization effect is absent (see subsection IV.D). Note 
that for the Au wire (E F = 5.6eV and m = 9.109 x 10~ 31 kg) we 
have A_p = 0.52nm, i.e., the wire width is W = 17.4Af = 9nm and 
the smallest A is only O.Olnm. Such small A is obviously not re- 
alistic, which demonstrates an important finding: The semiclassical 
Boltzmann result (full line) reproduces the classical scattering-matrix 
calculation (squares) for any Aa; only if A is too small to be realistic. 



produces the classical scattering-matrix calculation (squares) 
for any Ax only if A is too small to be realistic. This is 
another indication that the Boltzmann-equation theory of the 
edge-roughness-limited transport should be used with caution 
at least in the Au nanowires. 



V. SUMMARY AND CONCLUDING REMARKS 

We have studied quantum transport in quasi-one- 
dimensional wires made of a two-dimensional conductor 
of width W and length L 3> W. The main purpose of 
our work was to compare an impurity-free wire with rough 
edges with a smooth wire with impurity disorder. We have 
calculated the electron transmission through the wires by 
the scattering-matrix method, and we have obtained the 
Landauer conductance/resistance for a large ensemble of 



macroscopically identical disordered wires. 

We have first studied the impurity-free wire whose edges 
have a roughness correlation length comparable with the 
Fermi wave length. The mean resistance (p) and inverse mean 
conductance 1 / {<?) have been evaluated in dependence on L. 
For L — > we have found the quasi-ballistic dependence 
V(ff) = (P) = l / N c + PqbL/W, where l/N c is the funda- 
mental contact resistance and p q i, is the quasi-ballistic resis- 
tivity. For larger L we have found the crossover to the diffu- 
sive dependence l/(g) — (p) = \/N%** + pdifL/W, where 
Pdi f is the resistivity and 1 /N^f is the effective contact re- 
sistance corresponding to the N*** open channels. We have 
found the universal results p q b/ Pdif — 0.6N C and N*** ~ 6 
foriV c >l. 

Approaching the localization regime, we have demon- 
strated the following numerical finding: as L exceeds the 
localization length £, the resistance shows onset of localiza- 
tion while the conductance shows the diffusive dependence 
I /(g) ~ l/iV c e// + PdifL/W up to L ~ 2£ and the local- 
ization for L > 2£ only. On the contrary, for the impurity 
disorder we have found a standard diffusive behavior, namely 
l/(g) ~ (p) ~ 1/N C + PdifL/W for L < £. We have also 
seen that the impurity disorder and edge roughness show the 
universal conductance fluctuations of different size, as already 
reported in previous works HI 25 , 27 , 28 ] . 

We have then attempted to interpret our quantum-transport 
results in terms of an approximate but microscopic analyti- 
cal theory. In particular, the crossover from l/(g) = (p) = 
1/JVc + p qb L/W to l/(g) ~ (p) = l/N e Jf + p di fLlW in 
the wire with rough edges has been derived analytically as- 
suming the uncorrelated edge roughness and neglecting the 
localization effects. In spite of this approximation, our an- 
alytical results capture the main features of our numerical 
results. Specifically, we have derived the universal results 
p qb /Pdif = £N C and A^// ~ 2.5. 

We have also derived the wire conductivity from the semi- 
classical Boltzmann equation, and we have compare the semi- 
classical mean-free path with the mean-free path obtained 
from the quantum resistivity Pdif- For the impurity disorder 
we have found, that the semiclassical and quantum mean-free 
paths coincide, which is a standard result. However, for the 
edge roughness the semiclassical mean-free path strongly dif- 
fers from the quantum one, showing that the diffusive trans- 
port in the wire with rough edges is not semiclassical. We have 
found that it becomes semiclassical only if the roughness- 
correlation length is much larger than the Fermi wave length. 
In such case the resistance and conductance tend to scale with 
L /£ like in the wire with impurity disorder, also showing the 
universal conductance fluctuations of similar size. 

We end by a remark about our edge-roughness model. It 
is the same model as in the previous simulations 
2311 . We believe that most of the results obtained within the 



model are model-independent. For example, we have tested 
the correlation function T(q) of Gaussian shape and we have 
seen that the resulting transmissions remain quite similar to 
those in the figure ( fl~3b The only exception is the correlated 
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limit (1102l i. which predicts the dependence I oc Ax. This 
dependence is the artefact of our J-(q) choice. For example, 
for the Gaussian correlation function, the dependence I oc Ax 
is replaced by a much faster increase with Ax. On the other 
hand, the uncorrelated limit d84l > can be equally well derived 
for the Gaussian correlation function or for the exponential 
one, i.e., the result ( T84b is very general. 
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APPENDIX A: DERIVATION OF THE RELAXATION TIME 

At zero temperature the distribution d53l reads 
eFh, 



f n {k) = f(E nk ) 



-kr n (E nk )S(E nk - E F ) (104) 



and the Boltzmann equation d55l > can be written as 

^kd(E nk -E F ) = J2J2 w -y(^ k ') 

n' k' 

x [/«'(*/)-/„(*)]. (105) 
We set (1104b and d56| into ( 1105b . We obtain 



2tt, 



kS(E nk -e f ) = J2 EyK"' k '\ u \ n > fc)l 

n' k' 

S(E nk - E n , k ,)[k'T n ,(E n , k ,)6(E n , k , - E F ) 
kT n (E nk )5{E nk - E F ) 



(106) 



We proceed similarly to the Q2D theory [29]. First, we mul- 
tiply both sides of ( 1106b by k and sum over k. Second, on the 
left hand side we replace J2 k by (-^/2tt) J dk and integrate. 
Third, on the right hand side we replace T n (E nk )8(E nk —Ep) 
by T n (E F )S{E nk - E F ) and 8{E nk - E n , k ,)8(E nk - E F ) 
by 8(E n 'k' — Ep)5(E nk — Ep), which is justified due to the 
presence of 8(E nk — Ep). We obtain the equation 



mL 
^2 



n' k k' 

x S(E n , w - Ep) [k 2 r n {E F ) - kk'T nl {E F ) 
and from ( 1107b eventually 



(107) 



where K nn > is defined by equation 
the relaxation time d57l ). 



(108) 

From (1108b we obtain 



APPENDIX B: GOLDEN RULE FOR SCATTERING BY 
IMPURITY DISORDER 

If the scattering potential U is given by the impurity poten- 
tial 0, then 



JV 7 



\(n',k'\U\n, k)\l v = \(n' , k'\}^5(x-Xi)S(y-yi)\n, k) 



(109) 

Since the impurities are positioned at random, equation d 1 09b 
can be readily simplified as 



N, 



\(n', k'\U\n, k)\l = J2 \(n', k , \^5(x-x i )8{y-y i )\n, k)\\ Y 

(110) 

and the right hand side of dl 10b can be easy evaluated: 



|(n , ,fe'|Cr|n,fc)|2,= IjJV I 



1 Nl 



l 2 ni 
LW 



(HI) 

where nj — Nj /(WL) is the sheet impurity density. To ob- 
tain the right hand side of ( 111 lb . we have replaced the term 

WjE^iXl'iyijxKyi) by integral ^ J™ xl>{y)xl(v)dy. 
Setting (II 1 lb into (l58b we get the expression ( f60b . 



APPENDIX C: GOLDEN RULE FOR SCATTERING BY 
EDGE ROUGHNESS 

Assume first that the electrons are confined in the wire by 
the potential barriers of finite hight. Specifically, if the wire 
edges are smooth, the electron Hamiltonian reads 

ff0 = -^(^ + ^) + v_e(- y) + v + e(y-n 

(112) 

where V_ and V + is the hight of the potential barrier at y = 
and y = W, respectively, and 6 is the Heaviside step func- 
tion. If the wire edges are rough, the Hamiltonian reads 
H = H + U, where 

U = V_ [e(-d(x) -y)- Q(-y)] + 

+ V+[Q(y-h(x))-Q(y-W)} (113) 

is the perturbation potential due to the edge roughness, with 
d(x) and h(x) — W being the fluctuations of the edges (see 
figure[TJi. Now we also assume that |d(x)| <C W and \h(x) — 
W I <§C W and we approximate ( II 13b as 

U = -V-d(x)8{-y) - V+[h(x) - W]8(y - W). (114) 
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Then 



\(n',k'\U\n,k}\ 



V 2 

£,2 



/3=± 



dx / dye l(k - k ' )x 

U JO 



Xn' (y)Xn{y)df3{x)5fi(y) 



av 

(115) 



where<i + (a;) = d_(a;) = d(a;), <J+(y) = 8(y-W), 

and S-(y) = S(y). If we integrate in the equation ( II 15b over 
the variable y and then use the limit 



V- 



lim V_ X »(0) = hm V +X 2 n (W) = 



iW 3 



(116) 

both terms in the sum on the right hand side of ( II 15b can be 
rewritten into the form 



L nL 

d Xl / dx 2 e i{k - k 'X Xl - x ^{d ( Xl )dp(x 2 )). 

JQ 

(117) 



L 2 



where the brackets (...) symbolize the ensemble averag- 
ing (instead of the symbol av) and (dp(xx)dp(x2)) is the 
roughness-correlation function. The correlation function de- 
pends only on the distance \x\ — X2 1, so we rewrite it as 



(dp(xi)dp(x 2 )) = SpFpQxi - x 2 |), 



(118) 



where Fp(\x\ — x 2 \) is the normalized correlation function 
and 5+ and <5_ are the root-mean squares of the randomly 
fluctuating functions d + (x) and d-(x), respectively. We in- 
troduce the variable x 3 = x% — x 2 and simplify ( II 17b as 



-Oa 



L 2 



dxo 





— I— 5 ? 



L—x 2 



dx 3 e^ k - k '^Fp(\x 3 \) 



dx 3 e^ k - k '^Fp(\x 3 \), 



(119) 



To obtain the second line in dl 191 ), we have replaced the lim- 
its L — x-2 and — x 2 in the first line by oo and — oo, respec- 
tively. Except for very small x 2 , such replacement is justified 
because -P^da^l) decays with increasing x 3 to zero on the 
distance (correlation length) much smaller than L. 
The equation (II 15b can now be expressed as 



\(n',k'\U\n,k)\l v =J2 



A. n A.rif 



8 2 pFp{\k - fc'|), (120) 



where 



Mi) 



e i ^Fp(x 3 )dx 3 



(121) 



is the Fourier transform of Fp(x 3 ). Assuming the same ran- 
domness at both edges we can put 

6±=6, F ± (x 3 ) = F(x 3 ), F±{q)=T{q), (122) 



i.e., we skip the indices ±. We set ( 1122b into (1120b and we 
replace the symbol X)s=± i n ( 1 1 20t by the factor of 2. 

Finally, we set ( 11201 ) into d58l ). Performing summation over 
k and k' we obtain the matrix elements K nn > in the form d64b . 
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